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Abstract 


In this paper, we compute asymptotics for the determinant of the combinatorial Laplacian 
on a sequence of d-dimensional orthotope square lattices as the number of vertices in each 
dimension grows at the same rate. It is related to the number of spanning trees by the well- 
known matrix tree theorem. Asymptotics for 2 and 3 component rooted spanning forests 
in these graphs are also derived. Moreover, we express the number of spanning trees in a 
2-dimensional square lattice in terms of the one in a 2-dimensional discrete torus and also 
in the quartered Aztec diamond. As a consequence, we find an asymptotic expansion of the 
number of spanning trees in a subgraph of 1 ? with a triangular boundary. 

1 Introduction 

In this paper we study the asymptotic behaviour of the number of spanning trees in a discrete 
d-dimensional orthotope square lattice and in the quartered Aztec diamond. Let L(ni,... 
denote the d-dimensional orthotope square lattice defined by the cartesian product of the d path 
graphs Pm , i = 1,..., d, where nt, i = 1,..., d, are positive non-zero integers. We set rii = aiU, 
i = 1,..., d and write indifferently rii or throughout the paper. By rescaling the distance 
between two vertices on the lattice L(ni,... ,nd) with a factor of 1/n, the limiting object as n 
goes to infinity is a d-dimensional orthotope of size ui x • • • x ad, that we denote by Kd'- 


Kd ■■= [0,ai[x ■ ■ ■ X [0,ad[- 


The volume of Kd is 


d 




Let m S {1,..., d — 1}. An m-dimensional face of Kd is defined by 

{{xi,...,Xd) € Kd\ C {!,...,d} such that Xi^ G [0,ai,] , q=l,...,m 

and G {0, ai^}, q = m + 1,..., d}. 


The volume of the sum of all the m-dimensional faces of it'd is given by 


m 





l^ii<---<im^d q—1 
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For example, Vi is the perimeter and Vf the area of Kd- 

Asymptotics for the determinant of the combinatorial Laplacian on graphs have been widely 
studied, see for example [111 [3 min]. It is related to the number of spanning trees of a graph 
G, denoted by t(G), through the matrix tree theorem due to Kirchhoff (see [1]) 


r(G) 


1 


det^Ac 


where det*AG is the product of the non-zero eigenvalues of the Laplacian on G and |I^(G)| the 
number of vertices in G. In H, the authors developed a technique to compute the asymptotic 
behaviour of spectral determinants of the combinatorial Laplacian associated to a sequence of 
discrete tori. The technique consists in studying the asymptotic behaviour of the associated 
theta function which contains the spectral information of the graph. Consider a graph G with 
vertex set I^(G). For a function / defined on V{G), the combinatorial Laplacian is defined by 


^cfix) = - f{y)) 

y~x 

where the sum is over all vertices adjacent to x. Let {Xk}k denote the spectrum of the Laplacian 
on G. The associated theta function is defined by 

feev(G) 

To compute the asymptotic behaviour of spectral determinants on a sequence of d-orthotope 
square lattices, we express the associated theta function in terms of the theta function associated 
to the discrete torus with twice vertices at each side of the torus. This can be done because of the 
similarity of their spectrum. We then use the asymptotic results from |2] . The formula obtained 
relates the determinant of the Laplacian on the discrete lattice L(ni,... ,nd) to the regularized 
determinant of the Laplacian on the rescaled limiting object, which is the real d-dimensional 
orthotope Kd, and to the ones on the m-dimensional boundary faces of Kd, m = I,...,d — 1. 
Moreover, we compute asymptotic results for the number of rooted spanning forests with 2 and 
3 components. 

We will prove the following theorem. 

Theorem 1.1. Given positive integers exi, i = l,...,d, let det*A]^(^ain,...,adn) prod¬ 

uct of the non-zero eigenvalues of the Laplacian on the d-dimensional orthotope square lattice 
L{ain ,..., adu). Then as n ^ oo 


d—1 -| / «oo 

logdet*Ai(„^„_...,„^„) = CdV^n'^ - E ( / (1 “ e“'“) 

m=l ^ X-IO 


—At\d—m—2mtT \ Trd^ 

e hylt) -j j 


( 2 - 2 i-‘^)logn-H E E logdet*A, 

m—l 


oti-i X••• X ai. 


+ ^EC'r(-iriog(4m)+o(l) 


where det*AQi x---xai„ is the regularized determinant of the Laplacian on the m-orthotope aq x 
• • • X ai^ with Dirichlet boundary conditions. The constant Cd is 

Cd= r{e-* - e-^^*Io{2tf)^ 

Jo ^ 
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where /q is the modified I-Bessel function of order zero. 

Notice that the limiting object Kd can be decomposed in a disjoint union of m orthotopes, 
m G {0,1,...,d}. More precisely, let x ••• x {*q}^i C {l,...,(i}, denote the m- 
dimensional orthotope of side lengths 0 ;^^,..., which is open in R"*. Then 

d 

it:d = {0}u|J y ai, x---xaj„. 

m=l - 


This decomposition is reflected in the theorem by the appearance of the sum over this decom¬ 
position of the logarithm of the regularized determinant of the Laplacian on the m-dimensional 
faces, TO = 1,..., d. 

By expressing the eigenvalues of the Laplacian on the square lattice L{ni,n 2 ) in terms of the one 
on the two-dimensional discrete torus Z^/diag(2ni, 2n2)Z^, we derive a relation, which is stated 
below, between the number of spanning trees on these two lattices. 

Theorem 1.2. Given positive integers ni,n 2 , let T{L(ni,n 2 )) denote the number of spanning 
trees on the rectangular square lattice L(ni,n 2 ) and r(T(2ni, 2 ^ 2 )) the number of spanning trees 
on the discrete torus Z^/diag(2ni, 2n2)Z^. We have 


T(L(ni,n2)) 


_25/4r(T(2ni,2n2))^/^_ 

(nin2)i/4((3 -b 272 )’"! - (3 - 272)’^!) 1/2((3 + 2v/2)"2 - (3 - 2y2)"2)i/2 ' 


In [5], Kenyon computed asymptotics for spectral determinants on a simply-connected rec¬ 
tilinear region in R.^. Here we compute it in the particular case of a triangular region. More 
precisely, we consider the quartered Aztec diamond of order n, denoted by QADn, which is the 
subgraph of If with nearest neighbours connected induced by the vertices (fci,fc 2 ) such that 
^ 1+^2 ^ n and ki,k 2 ^ 0. Figure [T] illustrates QAD^. In [1], Ciucu derived a relation be¬ 
tween the characteristic polynomials of the rectangular square lattice and of the quartered Aztec 
diamond using combinatorial arguments. From this one can deduce a relation for the number 
of spanning trees. In the second part of this work, we present an alternative approach for it. 
Consequently, we derive the asymptotic behaviour of it, stated in the following theorem, which 
shows that it is related to the regularized determinant of the Laplacian on the triangle with 
Dirichlet boundary conditions. 


Figure 1: Quartered Aztec diamond of order 7. 


Theorem 1.3. LetT(QAD„) denote the number of spanning trees in the quartered Aztec diamond 
of order n. Then as n ^ 00 


log(T(QADj) 


2(^ Q 2“^ 

—- log (2 -b v^)n - - logn -b logdet*AA -b — log 2 -b o(l) 
TT 4 8 
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where G is the Catalan constant and det*AA is the regularized determinant of the Laplacian on 
the right-angled isoscele unit triangle with Dirichlet boundary conditions. 


1.1 Regularized determinant 


Let M be a Riemannian manifold with or without boundary and let Am be the Laplace-Beltrami 
operator associated to M. If M has a boundary we associate Dirichlet boundary conditions to 
Am- Denote by {Afcjfe the eigenvalues of Am- The spectral zeta function associated to M is 
defined for 3?(s) > dimM/2 by 


Cm{s) 


E 


1 



It admits a meromorphic continuation to the whole complex plane (see for example [13] for M 
with a boundary and [5] for the case M is a torus). The regularized determinant of Am can then 
be defined by 

logdet*AM = -Cm(0)- 


1.2 Preliminary result 

We prove the following lemma which will be useful in the next section to invert relations between 
the theta functions. Throughout this paper, we set an empty summation to be one by convention. 

Lemma 1.4. Let {iq\q^i he an increasing sequence of positive integers. Let {uqfq^i be a se¬ 
quence of positive integers and {ni^}q^i a subsequence of it. Let f,g : M. be two sequences 

of variadic functions such that for all I € N^i, 


fiui,,-■ ■ ,n„) E 9{na,,. ■ ■ ,nak). 


k=0 ai<"<ak 

Then the following inversion formula holds: for all I G 

i 

l-k 


k^O 


ai<---<ak 

K}S=iC{*,}Ui 


Proof. Let I G N^i. From relation o between / and g, we have 


( 1 ) 


E(-l)' 

k=0 

DdS=iCK}Ui 

"B-i)'-" y 


/c=0 




E E 


m—0 ai<---<art 


k—1 m—1 ji<---<jk ai<---<am 

Dd5=lC{i,}' 1 {adr=iCDdS=l 


( 2 ) 
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where in the second equality we used that 


fc=0 ji<---<jk k=0 

b,}S=iCK}Ui 

from the binomial theorem, where Cf = l\/{k\{l — fc)!). For I ^ k ^ m, the double summation 
can be rewritten in one summation as 

Z E = ctr i: . (3) 

{“<,}r=iCb9},=i 

Thus (HI) is equal to 

k—1 m—1 ji<---<jm 

{jq} q=l^{'^q} q=l 

I l—m 

m=l k=0 

{jq} q = l^{^g} q = l 

where the second equality comes from the fact that the only non-zero term in the summation 
over m is when m = I from the binomial theorem. □ 

Remark. If we set /(nq,..., ni,) = fi and g{ni-^ , ■ ■ ■, rii ^.) = g^ for all l,k € Njji, we recover 
the standard binomial inversion: 

i i 

fi=Y. ^^9k if and only if gi = A, for all I e N^i, 

k—0 k—0 

with fo = go = 1. 


2 Asymptotic number of spanning trees in the d-orthotope 
lattice 

2.1 Theta function 

The eigenvalues of the Laplacian on the square lattice L{ni ,..., Ud) are given by (see [5]) 

d 

{Afe}fc=o,i,...,Ar-i = {2d - 2Ecos(7rfci/ni), fc* = 0,1,...,- 1, for z = 1,...,d} 

i=l 

where N = nf=iThe d-dimensional discrete torus of size 2ni x ••• x 2nd is defined by 
the quotient Z'^/diag(2ni,..., 2nd)Z'^ with nearest neighbours connected. We denote it by 
T(2ni,..., 2nd)- The eigenvalues of the Laplacian on T{2ni ,..., 2nd) are given by (see [T^ l 

d 

{Afe }fe=o,i,..., 2 <iAr-i = {2d - 2 E cos(7rfcj/nj), ki = 0,1,... ,2ni - 1, for z = 1,..., d}. 

2=1 
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Notice that {A^}fc C {A^}fe. 

The theta function on the d-orthotope lattice L{ni, ..., n^) is given by 

"1 — 1 rid — l 

0L(„1,....",)(<)= E ■■■ E 

fci— 0 fed— 0 

and on the discete torus T{2ni,..., 2nd) by 

2ni —1 2nd —1 

0T(2„i....,2n.)W= E ■■■ E 

fci—0 Alci—0 

Therefore by expressing the theta function on the d-orthotope square lattice L(ni,... ,nd) in 
terms of the one on the torus T(2ni,..., 2nd), one can deduce the asymptotic behaviour from 
the results obtained in [5]. 

Let |a, 6] denote the set of successive integers {a, a + 1,..., 6}. In the theta function on 
L(ni,... ,nd), the summation is over the discrete d-orthotope |0,ni — 1] x • • • x |0,nd — 1] 
that we denote by Jd, while for the torus T(2n\,... ,2nd) it is over the discrete d-orthotope 
|0,2ni — 1| X • • • X |0,2nd — 11, denoted by Jd- We decompose Jd as a disjoint union of 
1-dimensional faces, I = 0,1,..., d. The 0 dimension is the point 0 € Z'^, we call it the 
root of Jd. For I G {l,...,d — 1}, the 1-dimensional faces are defined by a subset of Z'^, 
(fci,...,fed) C Jd, such that C {l,...,d} such that G Il,ni^ - 1], q = !,...,! 

and ki = 0, g = 1 -I- 1,..., d. We call the d-dimensional face the interior of Jd where no coor¬ 
dinate is zero, that is |1,ni — 1| x • • • x Jl,nd — 1]. For example, in the 2 dimensional case, J 2 
decomposes as: 

J 2 = {0} u (p, m - II X {0}) u ({0} X p, 712 - II) u (p, m - 11 X p, 772 - 11). 

In the theta function of the torus T{2ni, ..., 2nd), the summation is over 2'^ copies of Jd, namely: 

Jd = I_I { Jd + (ei, ■ • ■, ^d)) 

eiefo,";} 

where the unions are disjoint. Figure [2] illustrates the decompositions of Jd and Jd in the case 
d = 2. 
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(a) J2 




(b) J2 





Figure 2 : J2 and J2 with ni = 4 and 772 = 3 . The big black dot is the root of J2, the small black 
dots are the two 1 -dimensional faces and the white dots are the interior of J2. 
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Let define the theta-star function on the lattice L{ni -^,... , 71 ^,), with 1 ^ ii < • • • < i; ^ d, 
by the following expression 

Ui^—l 

E ■■■ E 

where the summation is over the interior of the Lorthotope of size x ■ ■ ■ x that is, the 
ki start at 1 instead of 0 for all i G {*i,... ,ii}. Therefore the theta function is related to the 
theta-star function by the relation 

d 

6'L(ni....,n,)W = E E (4) 

l—O 

To evaluate the theta function on T{2ni,... ,2nd) we sum over the roots of the 2^^ discrete 
orthotopes Jd, the ^-dimensional boundary faces of Jd and the interior of J^. Summing over 
the roots of the J^’s means that for all * = 1,..., d, ki is either 0 either Ui. There are = 
d!/(j!(d — j)\) number of ways to take j of the ki’s to be zero. In this case the corresponding 
exponential term of the theta function is The term / = 0 is then the sum of all the 

possibilities: 

3=0 

Then for each 1-dimensional boundary face of Jd, where I G {1,..., d— 1}, there are d—l of the fc^’s 
which are either 0 either ni, this positions the /-dimensional face in Jd. For this I, there are C^d-i 
number of ways where j of the fc^’s are zero and the exponential term is then And we 

sum over the interior of the /-dimensional face L(nij ,... ,ni^), where l^zi<---<b^d, which 
is by definition the theta-star function d2(rt n )(^) with a factor of 2* since this configuration 

appears 2* times. So for / G {1,..., d — 1}, we have 

j—0 l^ii<---<ii^d 

Finally we sum over the interior of Jd, that is, when all the fc^’s are greater or equal to one, which 
appears 2‘^ times. This gives the / = d term 

2 X(n,.....n.)W- 

The theta function on T{2ni,..., 2nd) is then the sum over the /-dimensional faces of the 2'^ 
orthotopes Jd'. 


^T(2«i,...,2nd)(0 - 2^^^ r -^- 'j ^ 

which is equivalent to 

(1 + e“'“)“'^^'T(2ni,....2nd)(0 = E^^^^ -1-6“"“)“' ^ ) (^) ■ 

Z=0 1^21 ^cZ 
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By setting 

/(m, ...,nd) = {l + e“'“)“'^ 6 »T( 2 ni..... 2 nd) (0 and g(n*i ,... ,ni,) = 2^(1 + 
in Lemma [Til it comes 

.„„,(() = ^ «TP„,.2„,„)(()^ 

m=0 jl<---<jm 

{jq} q=l^{^q}q=l 

From the above relation and relation a, the theta function on the d-orthotope square lattice is 
expressed in terms of the theta function on the d-dimensional torus: 

Oun .„,,(f) = E E y; . 2n.,-,(ty 

Rewriting the double multi-index summation in one summation using ([3]), we have 

/=0 m—0 l^ii<---<irri^d 

Therefore the theta functions are related by 

1 ^ 

= E ^T( 2 n.,...., 2 n.„)(i). (5) 

m—0 


2.2 Preliminary calculation 

Let {A^}j=o,...,Af-i be the eigenvalues of the combinatorial Laplacian on the square lattice 
L{ni ,..., Ud)- For small t > 0, the theta function on the torus 0T(2ni^,...,2ni^) behaves as 

m 

0T(2n.,,....2n.„)(i)-2™(n^2g)e-'™‘Lo(2t)™ = O(t), t^O. (6) 

9=1 

We follow the method derived in [2]. From relation a and the behaviour of the theta function 
at small t > 0 ©, we start by writing the theta function on the lattice L{ni, ..., n^) as 


E«-^^^‘ = E 


m—1 


1 — e 


-4t 


d—m 


1 — e” 


a /-I _ p-4t\ 

0 L(ng,...,n.)W- E e-2-‘Jo(2trKJn™-l 

m—1 ^ ^ 


to ensure the convergence of the integral of the Gauss transform that will appear below. By 
taking the Gauss transform of the above, that is, multiplying by 2se“® * and then integrating 


8 






with respect to t from zero to infinity, we have 


2s 




c2 _|_ 

j^O ^ ^ ^3 m-1 

^ d p 

s / 


1 — e 


-4t 


d—m 




“h 2s 


e-« t(i_e-4‘)d-m 


1 


E 






.)«) - - ^cr)* 


^oo d 

+ / e-"'‘((l - - 2^ + V CX{1 - e-‘^* f-"^)dt 

3o m=l 

For m € {1,..., d}, define the functions and ^ such that 

'1 — e~ 


(7) 


dsliis) = 2s e-^ ‘ f-^- ] e-2-‘Jo(2t)™dt 


( 8 ) 


and 


dsnijs) = 2s 

Jo 

E "tp,.. 

(9) 

Equation d?]) can then be written as 

Q ^ ^ -j ^OO 

E = E + E + ^2s / e-^^‘((2 - e-^r - 2‘^)dt. (10) 

m=l m=l -^0 

By integrating over s equations dH) and dS]) we get 

nUs) = {e-^ ‘(0T(2„,.....2n,)(t) - V,^{2nre-^^*Io{2tr - 1) + e"*)- 

and for m ^ d, 

^ p CO J± 

poo 

,Js) = -/ 

Jo 


dt 

t 


K 


/o 

X 


i: ...,.)«)-33^r>”6-“/„(2(r-Fcr)* 
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By integrating equation m above we get 

d d d 

^log(s2 + Aj^) = ^ V^n'^lt{s)+^ 'H^_„(s) + ^C^(-2)-''log(s2 + 4fc) + constant. (11) 

TYi—l m—1 k—1 

The asymptotic behaviour of the functions X'^ and „ as s —?► oo determine the constant of 
integration. We have 

Xf{s) = 2logs + o(l) and X^{s) = o(l) for m d as s —>■ oo 

and ^ 

= -^^logs + o(l) and H1 ^„(s) = o(l) for m ^ d as s oo 

and 

log(s^ + A^) =2{N — 1) logs + o(l) as s —>■ oo. 

1/0 

Therefore, equation m as s —>• oo yields 

2{N — 1) logs + o(l) = 27V logs — log s + (1 — 2^) logs + constant + o(l) 

so that the constant is zero. Evaluating equation m in s = 0 gives the logarithm of the product 
of the non-zero eigenvalues of the Laplacian on the d-dimensional square lattice T(ni,..., n^): 

d— 1 d d 

log ( n + E ^ nt,niO) + E (-2)'" log(4fc) (12) 

j^O m—1 m—1 k—1 

where 

Cd= r{e-^ - 

Jo ^ 

2.3 Asymptotic expansion 

By expanding (1 - in ^m,n(0)> it can be rewritten as 

.. d—m 

^l,n(0)=-^EClm(-lf 

k=0 

noo ni , 

^ E / ^^*(^T(2ni^,...,2ni„)(t) — ]^(2ai^ne ^‘/o(27)) — 1)-I-e *)-r^ 

for all TO = 1,..., d, where the e“* term is added to make the integral converge. It can be added 
since ^d-m(~^)^ = 0 for to = 1,..., d — 1. By splitting the sum over k we have 

d 

E = 

m—1 

1 _ roo ™ TO 

~^E E / (^r(2ni^,...,2ni„)(t) - ]^(2ai,ne ^‘/o(27)) - 1-|-e *) — 

m—1 l^ii<---<iyri^d ^ q—1 

1 ^ d-m „cjo rn , 

-^E E E / (e-"'=‘(0T(2n.,.....2n.^)(t)-n(2«*.^®”'‘^o(2t))-l) + e-‘)-. 

m=l k=l 9=1 
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From [5] Theorem 5.8], the asymptotic behaviour as n —>• oo of the k = 0 term is given by 

d „oo rn 

t 


1 u, „00 

~ ^ X/ X/ / (^T(2nij....,2rii^)(0 “ ^‘/q (2t)) - 1 + 6 ‘) 

m—1 q—1 

1 , ^ 

= H (21ogn-C4™/diag(2a,^.....2a,^)Z-(0))+o(l)- 


After a change of variable t n?t, the sum over the non-zero fc’s can be splitted as 

^ d—1 d—k 




k—1 m—1 


■ m j 

/ e~‘^^^^*{eT(2ni^,...,2ni^)(.n^t) - (2ai^(271^^))) — 

d 0 „_ 1 

'^r 

pOO 

1 ^—1 


(e 


—n'^t _ ^—Akn^t 


q^l 

’e-"'=”'‘(0T(2n.,....,2„.^)(n^i)-l)y 


dt 

T 


dt 


(13) 


From the propositions in [5J section 5], the first, third and fourth integrals tend to zero as n —>■ oo. 
The second integral tends to 


/°°(e-*-e-4'“'‘)^=log(4fc). 


The limit as n —>■ cx) of CSl) is then 
1 

'¥ 


- d—1 d—k ^ d—1 

-t^EE crcy„(-i)'-iog(4t) = -pi: (^fe(2d-fe_i)(_i)feiog(4^)_ 

k—1 m—1 




Therefore, the „ (0) term together with the constant term of equation (fT^ behave as 



E 

m—1 l^ii<---<im^d 


CR'^/diag( 2 o 


( 0 ) + ^EC',^-l)Mog(4fc) 

fc=i 


as n —>■ oo. 

We will now express the derivative of the spectral zeta function on the m-dimensional real torus 
R’"/diag(2 aij,..., 2ai^)Z™, that is CR™/diag(2ai 2 a )Z'"’ terms of the derivative of the 
spectral zeta function on the m-dimensional orthotope of size x • • • x . The eigenval¬ 
ues of the Laplace-Beltrami operator with Dirichlet boundary conditions on the m-dimensional 
orthotope of size x • • • x ai^ are given by 

m X ^ X 2 

Afc = E ( — ) ^ e (N*)'”- 

q=l 


So that the spectral zeta function on the m-dimensional orthotope of size ai-^ x • • • x with 
Dirichlet boundary conditions, denoted by Cctij x---xai„ j is given by 


Coii ! 




1 

7 T- 2 s 
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The spectral zeta function on the real torus ]R™/diag(2aij,..., is given by 


Ck™ /diag(2c(ij )Z'" (s) 


1 

(27r)2® 


E 



The spectral zeta functions are then related by 


iag(2Q;ij^ ,...,2QLi. 


‘(■s)— E^ E Ctiii X-'-Xa^j (s). 


i=l 

{jq}q = l 


Summing the above over all Zg, g = 1,..., m and over all m, m = 1,..., d, gives 


CK’^/diag(2cKij^ ,...,2ai^ )Z'^ (^) 

m—1 

^ d m 

= ^E E E^* E CaiiX-Xaj, (s) 

m—1 l^ii<---<im^d 1—1 ji<---<ji 

{3q}q = l^{'^q}^=l 

^ d m 

E Cttii x---xai, (s) 

m—1 1—1 l^ii<---<ii^d 


1 ^ 

4e e 


d 

E E 

m—1 l^ii<---<ir, 






where in the last equality we exchanged the sums over m and I and used the fact that 
— 2'^“! By expressing the derivative of the spectral zeta function evaluated in 
zero in terms of the regularized determinant of the Laplace-Beltrami operator on m-dimensional 
orthotopes, m = 1,..., d, with Dirichlet boundary conditions, we have 


1 

^E E 

m—1 l^ii<---<im^d 


iag(2aij^ ,...,2Q;i, 


-( 0 ) 


d 

E E logdet*Aa,^x-xai,„- 

m—1 l^i\<---<im'^d 


Putting everything together gives the asymptotic behaviour of the determinant of the Laplacian 
on the d-dimensional square lattice L{n \,..., Ud) 


-1 / ^oo 

logdet*Ai(„^„_...,„^„) = ^ f / (1 - 6 “'“) 

rr, = ^ \d0 


—4t\d—m^—2mt 


Wtrj ) 


d 

+ (2 - 2^“'^)logn+^ ^ logdet*A„,^x-xa. 

m—1 l^ii<---<im^d 

1 ^ 

+ ^E^™(-iri°g(4m) + o(l) 

m—1 


as n —>• oo. Notice that in the bulk limit the lead term is the same as in the case of the torus 
but lower order terms are deducted; for each m-dimensional face, m = l,...,d — 1, a term 
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proportional to V^n'^ is deducted. This can be explained by observing that in the torus case we 
have a periodic lattice while for the d-dimensional hypercubic lattice the periodicity is substituted 
by free boundary conditions. Spectral determinants of the limiting d-dimensional orthotope and 
each of its m-dimensional faces, m = l,...,d—1, also appear. Moreover the last constant term is 
new in the development. The terms in appearing from the boundary effect can be written 

in the following way: 


pOO 

/ 

^0 


— 4t\d—m^—2mt 


j d m „oo 7, 

Wtr- = C^d-mi-l)'^ / (e-‘ - 

* fe=o -^0 * 


These integrals are denoted by Jm{‘2k -I- m) in [5] by Glasser and are related to the Mahler 
measure of the hypercubic polynomial P{xi ,..., Xm) = 4:k + 2m + Y^^=i{xj by the relation 

m{P) = log 2 -I- Jm{‘2:k + to). In [B], Glasser calculates these integrals in terms of hypergeometric 
functions for to = 2 and 3. For to = 1, one explicity has 


4_d-l 




-4t\d-l-2t ] 


Io{2t)^ = ^y] Ctii-ir log(2fc + 1 + 2V^Tk) 


4^- 


d-l 

uE 


(see m Proposition 2.4]). 

Example. Consider the two dimensional rectangular grid ain x a 2 n. The volume of the limiting 
rectangle of size ai x a 2 is V 2 = aiQ :2 with perimeter Vi = 2q:i + 2 a2- From Theorem 11.11 it 
comes 

4G 1 3 

logdet*Ai(ai„,a 2 n) = —Vir? - -log{l + V2)V^n + -logn + logdet*xas 
-I- log det* Aqi + log det* Aaa - ^ log 2 + o(l) 
as n —> 00 , which is equivalent to the formula derived in [3 section 4.2]. 


2.4 Asymptotic number of rooted 2- and 3-spanning forests 

In this section, we derive asymptotics for 2 and 3 component rooted spanning forests in d- 
orthotope square lattices. Let denote the number of rooted fc-spanning forests on L{ni ,..., Ud) 
which is given by the fc-th power in of the characteristic polynomial: 


n 

3=0 



Following [3, by expanding the above in powers of {s/n)'^, one finds that N 2 and are related 
to Nf by 



The number or rooted spanning trees Af is related to the number of unrooted spanning trees 
T(T(ni,..., Ud)) by the relation 


d 

= (]4ni^T(L(ni,.. .,nd)). 
2 = 1 


13 






Recall that from equation m, we have that 


^log((s/n) 
j^o 


" + At) = 


d 

E 

m—l 


'H^_„(s/n) + Y Cd (“2) log((s/n)2 + 4fc) 


m—l 


k^l 


where 


We have 


For d'^i, 


Xi{s/n) = 


lim Xj{s/n) = Cd- 

n—^oo 


poo 1, 2 

hm n\xi{sln) - Cd) = lim nM (1 - = ^-Wd 

n—foo n—^oo Jq t Z 

where Wd is the so-called Watson integral for the d-dimensional hypercubic lattice 

nOO 

Wd= / 

Jo 

In [ 8 ], Joyce and Zucker introduced the generalised lattice Green function defined by 

1 

Gd{n;k,w) = -^ TT(t)dt 

i Jo 

where n = {ni,... ,nd} is a set of non-negative integers, w ^ d, k > 0 and F is the gamma 
function. Here the lattice Green function will only appear with n = 0, hence we denote it shortly 
by Gd{k,w). In [S], numerical evaluations of the integrals Cd and Wd are computed and also in 
for Cd and in [7] for Wd and GdO,w). 

For d^ 5, 

lim {n\X^{s/n) - Cd) - n^s^Wd) = lim T (l- - 4d e-^‘^*Io{2t)‘^- 

n^ca n—>oo Jo \ 71^ J t 

= -^G,(2,d). 


Continuing in this way, we arrive at the following expansion for Xf^s/n) as n —>■ oo 


L(d-1)/2J 


k=l 


2k 

n‘^Xj{sln) = Cdu'^ + Y d{k, d) + o(n). 


(14) 


Recall that for m S {1,..., d — 1}, 

-I p CO J± 


We have 


hm Xi{s/n)=X^{0). 

n—^oo 


14 




Similarly as for m = d, we obtain as n —> oo 




L(m-1)/2J 


+ 


/^d-r, 


E (-1) 


k^m — 2k 


k^l 


'2k poo 

^ (1 - + o{n). 


(15) 


The above integral can be expressed in terms of the generalised lattice Green function: 

-| noo d—m 

/ (1 _ ^ ^ C'_„(-l)'G^(fc, m + 21). 

{k-l)Uo 

Putting equations (fHl) and (fT51) together gives the expansion for d ^ 3 

d-l L(d-1)/2J 2fc / 

Y,^og{{s/n)^ + Xf)=V,^c,n^+Y,V:ili{0)n^+ ^ _ / (_l)'=+iv(fG,(fc, 

j^O m—1 k—1 ^ 

/_1 d—1 

I j: \ ^) \ ^ -[rd^m — 2k 

+ IP 

' m=2fc+l 

where Sd^do = 1 if d ^ do and 0 otherwise. 

For d = 3, 

^ log((s/n )2 + = Vic3n‘^ + V,^l!{0)n + Vill{0)n^ + ^ViW^n + o{n) 

j/o 

with the special values [Wz is given in [S]) 

Ii'( 0 ) = ^log((17 + 2V2)(5-2y6)) and VF 3 = ^(^3 - 1 ) (^(l/24)^(ll/24))^ 
For d ^ 4, 


/^d-r 


( 1 -e-n 


i—m^—mt ] 


^dt ) + o{n) 


d-l 


J2 log((s/n)2 + A,^) = Vd^Cdn‘^ + ^ P^T;l(0)n" 




2 / d,-l .. .00 \ 

+ y \VfWdn‘^-^ - J2 I (1 - 


4 / d-l .. .00 \ 

+ y \-Sd^zVd^Gd{2, d)n'^-4 + E ^(1 " e-^‘)"-’"e-™*/o(t)'"dt j 


+ (terms in with fc ^ 3) + o{n) 


as n —?► 00 . 


On the other hand, the formal expansion of log(('S/n)^ + Xj) gives 

(-l)P-i 1 


0^0 


Y. iog((./n)^+Ap=log (n Af)+E ^ (;) i: 




p^l 


’ 
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By identification of the terms in we find the asymptotic number of rooted 2-spanning forests 
for d = 3, as n —> oo 

Ni = + o(n^)^ 

and for d ^ 4, we have 


1 wd ^_1 Y^- 1 r°° 

J2jl= - E / (1 - ‘/o(t)'"dt + o(n3) 

iVO J ™=3 


SO that 


C yd d, 1 . -oo \ 

-fWdn‘^ - ^ (1 - e-^Y-^e-^^Io{trdt + oin^)j 


where Nf is asymptotically given by Theorem ll.il By identification of the terms in s^, we find 
that as n —?► oo 


E' 


1 


..0 


= O(n^) 


so that the asymptotic number of rooted 3-spanning forests for d = 4 is given by 


Nl = 


T/ 4 T /4 

.^8 _ 


POO 

/ 

do 


e ‘)e ‘^*Io{tydt + o{n")]N{, as n —>■ oo. 


Remark. It would be interesting to find the next terms in the development. For the 2-dimensional 
case, we would need to find the asymptotic development as n —>■ oo of the following integral 

„2(i _ e-("/'^)'‘)e-^*/o(2t)2*. 



In [3, the authors computed the asymptotic development of '^j^Q^og{{s/n)‘^ -(- Xj) in the case 
of the torus with other techniques. To generalise their result to higher dimensions with our tech¬ 
niques, for example in the 3-dimensional case, we would need to find the asymptotic development 
of ^ 

/'°°(n3(l - - nsH)e-^^h{2tf — 

do ^ 

as n —>• oo. This would enable us to derive asymptotics for the number of rooted A:-spanning 
forests with fc > 4. 


2.5 Spanning trees in two-dimensional square lattices 

In the two-dimensional case, one can derive an exact relation between the number of spanning 
trees on the rectangular square lattice ni x n 2 and the one on the torus of size 2ni x 2n2. The 
product of the non-zero eigenvalues on T(2ni, 2n2) is given by 


2 ni —1 2 n 2 — 1 

det*Aj’(2„j_2n2) = n n (4 — 2 cos(7rfci/ni) — 2 cos(7rfc2/n2)). 

fcl —0 ^ 2—0 

{ki,k2)^0 
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The product over fci, k 2 is a disjoint union of products over four squares of size ni x 712 - We split 
this product as a product over the 0, 1 and 2 dimensional faces of the squares. It comes 

n\ — l 722 — 1 

det*AT( 2 ni, 2 n 2 ) = 4^8 (2 - 2cos(7rfci/ni))^ (2 - 2 003 ( 77 ^ 2 / 712 ))^ 

ki^l k2^l 

721 —1 722 — 1 

X ^ (6 — 2 cos(7rA:i/ni))^ (6 — 2 cos(7rA:2/n2))^ 
fci^l fc2^1 

721 —1 72 2-1 

X n n (4 — 2 cos(7rfci/7ii) — 2 008 ( 77 ^ 2 / 712 ))'^. (16) 

fel = l /C2 = l 

On the other hand, the product of the non-zero eigenvalues on the square lattice L(ni,n 2 ) is 
given by 

721 — 1 722 — 1 

det*Ai(„j^„2) = II ^ 003 ( 77 / 71 / 771 ) - 2 003 ( 77 / 72 / 712 )). 

All—0 A:2—0 
(fcl.fc2)#0 

By splitting the above product as a product when ki = 0, then k 2 = 0, then 1 ^ /71 ^ 711 — 1, 
1 ^ /72 ^ 712 — 1 , we get 

721-1 722 — 1 

det*AL(„i,„ 2 ) = (2 - 2cos(7r/7i/7ii)) ( 2 - 2003 ( 77 / 72 / 712 )) 

/ci —1 A;2 —1 

721 — 1 722 — 1 

X n n (4 — 2 003 ( 77 / 71 / 711 ) — 2 003 ( 77 / 72 / 712 )). (17) 

All—1 A:2 —1 

Using the matrix tree theorem and putting the following identities coming from relations for 
Chebyshev polynomials of the second kind 


n (2 - 2 cos(^fc/77)) = 77 and [f (6 - 2 003 ( 77 ^/ 77 )) = (^ + ^^)" ~^3 - 2a/2)» 

Ai^l k^l 

in equations (fT51) and m, it follows that 


T(L(r7i,r72)) 


_25/V(r(277i,2T72))^/^_ 

(77i772)i/ 4((3 -f 2^2)"! - (3 - 2y2)"i)i/2((3 + 2^2)"= - (3 - 2V2)^^y/^ ' 


Remark. It would be interesting to see if one could generalise the above relation to higher 
dimensions. It could not be done in the same way as it is done above. More precisely, when 
splitting the product in 3 dimensions as a product over 0, 1, 2 and 3 dimensional faces, one would 
need for example to evaluate the following product 


721-1 722 — 1 

n n (8 — 2 cos(7r/7i/r7i) — 2 cos{'Kk2 / n2)) 

All —1 A:2 —1 


appearing for the 2-dimensional face defined by k^ = 773 and ki = 1,..., 771 — 1, k 2 = 1,..., 772 — 1. 
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3 Asymptotic number of spanning trees in the quartered 
Aztec diamond 


3.1 A relation between the number of spanning trees on the quartered 
Aztec diamond and on the square lattice 

In [H IIP) , the authors showed that the number of spanning trees in the quartered Aztec diamond 
of side length n is given by 

r(QADj = n (4 — 2 cos(7rfci/n) — 2 cos{'Kk 2 / n)). 

{)<k\<k2<n 

The product of the non-zero eigenvalues on the square grid of side n is given by 

n—1 n—1 

det*AL(„,„) = n n (4 — 2 cos(7r/ci/n) — 2 cos{'Kk 2 / n)). 

ki —0 k2 —0 
(fcl.fc2)#0 

By splitting this product as a product when fci = 0, then ^2 = 0, then ki = k 2 , ki = 1,..., n — 1, 
then ki < k 2 and ^2 < fci, we have 

n—1 n—1 

det* (2 — 2 cos(7rfc/n))^ (4 ~ 4 cos^nk/n)) 

k=l k=l 

X (4 — 2cos(7rfci/n) — 2 cos(7rfc2/n))^. 


From the matrix tree theorem, it follows that 


t(QADJ 


T{L{n, 
Vii2("-i)/2 ■ 


(18) 


3.2 Asymptotic expansion 

From (la, we have 


log{T{QADn)) = ilogdet*Ai(„_„) - ^ log2 - ^ logn-b ^ log2 (19) 

where the asymptotic behaviour of logdet^A^^^ is given by 

4:0 3 1 

logdet*Ai(„_„) = —- 21og(l + ^2)71-b - logn - C(xi( 0 ) “ 2 C)( 0 ) - ^log2 -bo(l) (20) 

as n —>■ 00 . 

Consider the right-angled isoscele triangle with the sides of same length equal to 1. The eigen¬ 
values of the Laplace-Beltrami operator with Dirichlet boundary conditions are given by 

Afe = '^'^{kl + k^), with k = (fci, ^ 2 ) € (N*)^ and fci > k 2 . 

The associated spectral zeta function with Dirichlet boundary conditions, denoted by Caj is then 
given by 


Ca(s) 


— Y _ 
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The spectral zeta function on the unit square with Dirichlet boundary conditions is given by 


Cixi(s)-^ E - 2 Ca(s )+2 XlW. 

The spectral zeta function on the unit interval with Dirichlet boundary conditions is related to 
the Riemann zeta function by Ci('S) = (2/7J‘^®)C(2s) with special values in 0, C(0) = —1/2 and 
C'(0) = — (1/2) log(27r). Thus we have that Ci(0) = —2 log 2. By differentiating the above and 
evaluating in s = 0, we get 

C;xi(0) = 2Ck(0)-log2. (21) 

Putting (IT^ . (I^Tll together and writing the derivative of the spectral zeta function in 0 

in terms of the regularized determinant of the Laplace-Beltrami operator on the right-angled 
isoscele unit triangle with Dirichlet boundary conditions, that is 

Ca( 0) = -logdeDAA, 

gives the asymptotic behaviour of the number of spanning trees on the quartered Aztec diamond, 
namely 

2C 3 21 

log(T(QAD„)) = —- log(2 -|- V2)n - - logn -|- logdet*AA -f — log2 -|- o(l) 

TT 4 O 

as n —?► oo. 
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